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This letter discusses a novel use of the Geant4 simulation toolkit to model molecular transport in 
a vacuum environment, in the molecular flow regime. The Geant4 toolkit was originally developed 
by the high energy physics community to simulate the interactions of elementary particles within 
complex detector systems. Here its capabilities are utilized to model molecular vacuum transport in 
geometries where other techniques are impractical. The techniques are verified with an application 
representing a simple vacuum geometry that has been studied previously both analytically and 
by basic Monte Carlo simulation. We discuss the use of an application with a very complicated 
geometry, that of the Large Synoptic Survey Telescope camera cryostat, to determine probabilities 
of transport of contaminant molecules to optical surfaces where control of contamination is crucial. 
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I. INTRODUCTION 

The cryostat of the Large Synoptic Survey Telescope 
camera (LSSTCain}llJ represents a significant challenge 
for vacuum design 0. With a volume of 2.9 m'^, it will 
contain nearly 1000 kg of material, and have more than 
40 m^ of exposed surface area within. The geometry 
of the LSSTCam cryostat presents an environment too 
complicated for vacuum transport to be modeled by ana- 
lytical calculations, or by techniques such as thermal con- 
ductance calculations in conjunction with finite element 
analysis [sj. Instead, a fully three dimensional model of 
the geometry is required, with a full Monte Carlo treat- 
ment of the initial conditions and subsequent propagation 
of particles through their interactions with surfaces. Such 
calculations have been performed previously for relatively 
simpler geometries^. The complexity of the LSSTCam 
cryostat would render the specification of the geometry 
by similar means impossible. Therefore, the Geantd sim- 
ulation toolkit [Bl, ^ presents itself as an ideal platform 
for developing an application to model in fine detail vac- 
uum transport in the LSSTCam cryostat and similarly 
complex vacuum structures. 

Geantd is a Monte Carlo simulation toolkit used ex- 
tensively in particle physics, nuclear physics, and medical 
physics applications. It is a flexible, self-contained soft- 
ware package that allows the definition of a volume con- 
taining real materials arranged into complex geometries 
and then the transport of elementary or composite parti- 
cles within the volume. Geant4 follows any fundamental 
interactions of the particles with the materials, as well as 
any products of those interactions. The toolkit consists 
of a large number of C++ class bases which the applica- 
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FIG. 1: CAD rendering of an LSST camera cryostat design in 
cross section (this particular scenario is known as "option 3"). 
A brief summary of the cryostat is given in 21and its geometry 
is discussed in ijlVI Shown here is one science raft tower mod- 
ule containing CCD detectors, electronics boards and other 
structures, the grid (rendered in blue), the cryoplate (pink 
and dark orange), the cold shroud (red), the lower shroud (or- 
ange) , some of the pumping plenum and chimneys (dark green 
and white), the feed-thru flange (transparent light green), and 
the pump plate (blue). The cryostat body and the L3 lens 
are rendered in grey. Some structures are not shown. 



tion developer uses in order to create a specific applica- 
tion that models a physical system. To develop a Geantd 
application, particles, their interactions, a detector ge- 
ometry, and relevant cataloging of outcomes (known as 
"scoring") must be implemented. 

In this letter we describe Geantd applications that were 
developed to model vacuum transport in the molecular 
flow regime. The general strategy is summarized in 
In mill an application with a simple geometry to verify 
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the techniques used is discussed. In mV\ a specific ap- 
pUcation to the LSSTCam cryostat is described. Several 
cases of the cryostat geometry and initial conditions of 
the molecular distribution are evaluated to demonstrate 
the power and flexibility of this design tool. These gen- 
eral techniques could of course be applied to any vacuum 
geometry, and this is explored in ifVl^ 



II. APPLICATION STRATEGY 

To model molecular transport in vacuum with the 
Geant4 toolkit, the vacuum space of interest is de- 
fined as the detector geometry. All structures are spec- 
ified according to their material (stainless steel, cop- 
per, GIO, silicon, and so on) in the standard Geant4 
way taking advantage of the material database. They 
are also assigned their predicted operating temperatures. 
Custom particles are specified with implementations of 
G4ParticleDefinition, to represent the desired molecu- 
lar species ("Waterons" and "Nitroginos" for example, 
representing water and nitrogen molecules, respectively) 
that lack internal structure but have the proper molecu- 
lar weight, allowing, for example, the mean velocities and 
residence times to be affected by contact with surfaces of 
different temperatures. 

TABLE I: Test case described in ijllll — probability of a parti- 
cle first reaching the opposite end rather than the originating 
end of a pipe of length L and radius R via molecular vacuum 
transport. is the probability as determined by analytical 
calculations and reported in Clausing[3l, Pm is the probabil- 
ity as determined from Monte Carlo methods by Davis 4], and 
Pa is the probability as determined by a Geant4 application 
described in ijllll with 10000 randomized runs. 

L/R Pe Pm Pg' 

0.5 0.80 0.80 0.80 

1.0 0.67 0.68 0.67 

1.5 0.58 0.57 0.57 

2.0 0.51 0.53 0.51 

3.0 0.42 0.43 0.44 

4.0 0.36 0.36 0.35 

5.0 0.32 0.32 0.29 

10.0 0.20 — 0.18 

20.0 0.11 — 0.10 

There are no particle-particle interactions, as the phys- 
ical situations being modeled are high vacuum systems 
where transport is in the molecular flow regime, so that 
the "ballistic" approximation is valid. For example, the 
LSSTCam cryostat will have pressures of ~10~^ torr 
upon pump-down and around two orders of magnitude 
lower when surfaces are cooled to their operating tem- 



AII applications discussed in this work are available at 
|htt p : / / www . slac . Stanford, edu / ~j acks [ 



peratures, so that the molecular mean free path is sig- 
nificantly larger than the dimensions involved. Particles 
thus only interact with material surfaces, and then in 
only a particular way that represents how surfaces direct 
particle transport in the molecular flow regime. When 
a particle encounters a surface, it is adsorbed at the 
point of contact and then desorbed with the probabil- 
ity of a given direction proportional to the cosine of the 
angle with the normal vector to the surface at the point. 
We specify these interactions in an implementation of 
G4VDiscreteProcess. 

An implementation of G4GeneralParticleSource is 
used to originate particles with initial momentum direc- 
tions and positions chosen randomly according to dis- 
tributions appropriate to the application, and from the 
regions of interest within the volume of the vacuum. 
Particles are tagged with the name of the first sur- 
face they encounter, which can be considered the orig- 
inating item, and this is done in an implementation of 
G4UserSteppingAction with the originating surface name 
and any other desired information attached with an im- 
plementation of G4UserTrackInformation. 

Scoring is implemented such that particles propagate 
until they reach one of several particular specified fates, 
for example exiting through a pump port, or contacting 
a particular surface. These terminating fates are im- 
plemented through the use of regions specified in the 
detector construction, which are checked at each step 
in our implementation of G4UserSteppingAction. Each 
molecule's originating item and terminating fate is cata- 
loged in the output data file, which consists of columns 
representing the different fates with a "1" entered into 
the proper fate column for each molecule. In this man- 
ner the total number of molecules resulting in each fate 
can be easily summed. 



III. CONFIRMATION OF TECHNIQUES WITH 
A SIMPLE GEOMETRY 

To test the validity of using the Geant4 toolkit to 
model molecular transport in vacuum with the tech- 
niques discussed here, we have developed a Geant4 appli- 
cation with a very simple geometry for which transport 
results can be compared to those in the literature. Fol- 
lowing Davis [3] , the geometry modeled is a right circular 
cylinder, where molecules enter at one end with randomly 
distributed initial radial and azimuthal positions in the 
plane of the opening and momentum vectors randomly 
distributed proportional to the cosine of their angle with 
the surface normal of the opening plane. Molecules are 
propagated through the cylinder with wall interactions 
as discussed in ^TT] The quantity of interest is whether a 
molecule first reaches the opposite end of the cylinder or 
returns to the opening through which it entered. Thus 
this is just the problem of molecular transport probabil- 
ity through a pipe of a given length and radius. 

The probabilities vary with the ratio of the length 
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FIG. 2: View of an OpenGL rendering of the "option 3" ge- 
ometry used in the LSST camera cryostat Geant4 application. 
For clarity, certain structures are rendered in wireframe and 
others are rendered invisible. 

(i) to the radius {R) of the cylinder, and have been 
calculated analytically by Clausing^?] in 1930 and sub- 
sequently by Davis [4] using simple Monte Carlo meth- 
ods. A Geant4 application for this geometry is imple- 
mented and the results for various L/R values compared. 
The application has regions specified at the two ends of 
the cylinder that terminate a given molecule's trajectory 
when they are entered, and the terminating region is 
recorded in the output file as described in ^TT] Table 1 
shows the probabilities for a molecule to reach the op- 
posite opening for different values of L/R, as calculated 
analytically by Clausing, via Monte Carlo methods by 
Davis, and by the Geant4 application. The statistical 
uncertainty on the probabilities we have obtained can 
be estimated from standard theory as tr = [y/n/M]/M, 
for a situation where there are n randomized trials with 
M possible discrete outcomes. In this case, cr=0.005. 
The results from the Geant4 application show very good 
agreement with the probabilities obtained by the other 
methods. 



IV. LSSTCAM CRYOSTAT APPLICATION 

Figure [1] shows a CAD rendering of an LSSTCam cryo- 
stat layout. The strict requirements on LSST's photon 
throughput necessitate a detailed knowledge of the effect 
of potential molecular contaminants that may deposit on 
the cold CCD focal plane surface within the cryostat. 
One crucial input to any model that addresses this ques- 
tion is the transport of molecules within the cryostat. In 
particular, how likely is a contaminant molecule from the 
bulk of the vacuum space where most of the electronics 
and other structures are located to reach the focal plane 



region rather than being removed from the system else- 
where? 

TABLE II: Likelihood of three different outcomes (exiting via 
a pump port, encountering a getter pump, or encountering the 
front surface of a CCD) occurring first for molecules in the 
LSSTCam cryostat application, for several different geome- 
tries. The geometries are discussed in ijlVI Results are shown 
both for molecules originating with random positions and mo- 
mentum directions from the entire region within the cryostat 
behind the focal plane ("All loc"), and for molecules origi- 
nating with random positions and momentum directions only 
within the envelope containing all of the electronics boards 
("Elec. reg."). Values may not total 100% due to rounding. 
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24% 
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50% 


17% 


33% 



A brief description of the LSSTCam cryostat geometry 
follows: The cryostat body, which forms the vacuum en- 
closure, is a truncated conical section with small interior 
dimension (ID) 0.95 m and large ID 1.06 m. It is sealed 
at the smaller "front" end by a glass lens designated 
"L3". Directly behind L3 is the focal plane consisting 
of the 3.2 billion 10 CCD pixels arranged in 21 raft 
tower modules (RTMs). Each RTM contains nine 4kx4k 
pixel CCD detectors, and is itself a self-contained cam- 
era with the CCD packages, a Silicon Carbide "raft" on 
which they arc mounted and aligned, conductance barri- 
ers, electronics boards, thermal and wall structures, and 
connectors and cabling. Four triangular shaped RTMs 
carrying guide and wavefront sensors are located at the 
corners of the focal plane. The RTMs pass through a Sil- 
icon Carbide ceramic honeycomb structure known as the 
grid, which kinematically supports the raft. Behind the 
grid, the cryoplate is a mechanical structure that carries 
the balance of the load of the RTM. The cryoplate con- 
tains cryogenic refrigerant channels that provide cooling 
to the RTM electronics and other structures. Behind the 
cryoplate, a separate cold circuit removes heat from the 
rearmost electronics at a higher temperature. There are 
also various shrouds, chimneys, and plenums that provide 
thermal radiation shields and direct molecular transport. 
The rear of the cryostat is an annular feed-through flange 
containing hermetic signal and cryogenic feed-throughs, 
and an octagonal plate containing ports for turbomolecu- 
lar and ion pumps and gauges. In combination, they seal 
the rear of the vacuum space. Two large molecular sieve 
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FIG. 3: View of an OpenGL rendering of the "baseline" geom- 
etry used in the LSST camera cryostat Geant4 apphcation, 
along with an example trajectory of a simulated molecule. 
This particular molecule has undergone a far lower than av- 
erage number of surface interactions and ended up contacting 
the front surface of a CCD. 



getter pump structures reside just behind L3, facing the 
focal plane. Additional large area cold activated-charcoal 
getter pumps will be placed in the transport regions be- 
tween the front and the rear of the cryostat, but are not 
yet included in the current models. 

Three different primary cryostat layouts have been un- 
der consideration for LSSTCam. The first, or "baseline," 
design features RTMs with distinct front and rear elec- 
tronics boards connected by 2m long cables. Inside the 
RTM, the frontmost boards are interspersed with copper 
thermal planes and connected to the cryoplate. An addi- 
tional warmer coldplate structure is provided to support 
and cool the rearmost boards. A modified design, "Op- 
tion 2" contains the electronics on a single longer board 
with built-in copper thermal bars, replaces the coldplate 
with a simpler tubular structure, and in addition makes 
several other modications. "Option 3" modifies the de- 
sign even further by reconfiguring the radiation shrouds 
and pumping plenum geometry. To explore the design 
parameters, further modifications to Option 3 have been 
examined, for example increasing the diameter of the ion 
pump ports by either two or four times. These cases 
are designated "Option 3 -I- 2x1" and "Option 3 -I- 4x1" 
respectively. Figures [2] and [3] show examples of the detec- 
tor geometry developed for this application as rendered 
in the OpenGL viewer. 

For this application scoring has been implemented such 
that each particle interacts with surfaces until it reaches 
one of three fates — either contacting the front surface 
of the CCDs in the focal plane, contacting the surface of 
a getter pump located near the focal plane, or exiting the 
cryostat through one of two ion pump ports. These three 



fates, any of which ends any given molecule's trajectory, 
are determined through the use of regions specified in the 
detector construction as discussed in ^|TT1 As an indica- 
tion of the complexity of the LSSTCam geometry, it is 
found that molecules can have anywhere from less than 
one hundred to hundreds of thousands of surface encoun- 
ters before encountering one of these three terminating 
regions. 

The probabilities for each of these outcomes for several 
alternate geometries under consideration are presented in 
Table II; each is based on simulations of 1000 particles. 
The statistical uncertainties on these values can be ob- 
tained as outlined in ijllll and are a=l%. These results 
show that in the first three geometries most molecules 
(>85%) eventually reach the focal plane region rather 
than exiting through the rear ion pump ports. Once 
there, they then encounter either the getter pump or the 
front surface of the CCD detectors. From this one con- 
cludes that the getter pumps will need to have a large 
enough capacity to handle this fraction of the flow. It 
also suggests these pump surfaces should be made larger 
or be better positioned. Alternately, modifications to 
the geometry to make molecular transport into the ion 
pumps more likely is warranted. The simplest modifica- 
tion would be to increase the diameter of the ion pump 
ports, and as shown in Table II this by itself greatly in- 
creases the likelihood of a molecule exiting through an 
ion pump port before reaching the focal plane region. 
The ability to readily test the effects of such modifica- 
tions makes a strong case for this Geant4 modeling and 
simulation approach. 

A number of features within the cryostat have not yet 
been rendered for the application detector geometry, be- 
cause they would be too cumbersome to implement or 
are not yet specified, nor are they likely have a signifi- 
cant effect on the bulk molecular transport issues under 
investigation. These features include fasteners, cables, 
thermal straps, fluid lines, hermetic feed-throughs, CCD 
package kinematic mounts, raft hold downs, and perime- 
ter cutouts in the grid. Additionally some features have 
necessarily been simplified to allow for their specification 
with Geant4 detector construction. It is unlikely that 
these simplifications and omissions would have a signi- 
cant effect on the bulk vacuum molecular transport issues 
being modeled. 



V. DISCUSSION 

Geant4 provides a flexible and well-documented frame- 
work to simulate particle transport and interactions 
across complex geometries. These features, in concert 
with the techniques described herein, provide a unique 
method to model molecular transport in complex vac- 
uum environments, where analytical calculations or other 
techniques are impractical. Such is the case in many 
real applications like the LSSTCam cryostat. As demon- 
strated in mm these general techniques reproduce results 
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FIG. 4: View of an OpenGL rendering of the geometry of a 
Geant4 application representing a section of accelerator pipe, 
where regularly spaced irises form cavities, as discussed in ^V] 
A trajectory of a molecule created with initial position on the 
wall of the middle section and with random initial momentum 
is shown. 

from analytical calculations and other Monte Carlo meth- 
ods that are available for comparison. 

The techniques we discuss are clearly not unique to a 
particular vacuum geometry and could be used to model 
molecular transport in any geometry. Indeed such mod- 
eling can also provide additional information as well. As 
an example. Figure S] shows a rendering of the geom- 
etry of another Geant4 application we have developed, 
that of an electron accelerator beam pipe with regularly 
spaced irises that form cavities, with pumps located at 
either end. This application has been used to determine 
mean residence lengths of molecules in a given cavity sec- 
tion, and therefore the partial pressures within, given an 
outgassing rate for the interior surfaces of the pipe. To 
sum the path lengths within a given volume, one simply 
carries the net length through in an implementation of 
G4UserTrackInformation and adds to it in an implemen- 
tation of G4UserSteppingAction, where the step length 
is added if the step terminates within the volume of in- 
terest. 

As discussed in mVl for the LSSTCam cryostat appli- 
cation, we have so far implemented a simple scoring (i.e.: 
fraction of particles reaching a focal plane structure ver- 



sus reaching a pump first). As already indicated, the ef- 
ficacy of simple changes to the model such as increasing 
the diameter of the ion pump ports can be readily eval- 
uated. In further studies, the effects of adding, resizing, 
and moving pump ports and getters will be evaluated. 
More extensive changes to the geometry and design (e.g.: 
plenum shape, charcoal getters, molecular flow barriers, 
etc.) will also be studied. 

Thus far, interest has been in the steady state or 
"worst-case" scenario estimates. However these appli- 
cations can be easily extended to include temperature- 
dependent vacuum sticking coefficients representing the 
mean duration for a molecular species to be adsorbed 
on a particular surface. This time dependency can be 
included in evaluating, for example, whether a particle 
will actually contribute to building monolayers of con- 
taminants on the front surface of the CCDs within a 
given interval of time. Results from this application are 
crucial for both evaluating the transport, and therefore 
focal plane contamination consequences, of different lay- 
outs and features under consideration, and for providing 
a necessary component of the mapping from outgassing 
of components (that are measured in a test stand 0) to 
the resulting effect on the contaminant levels adsorbed 
on the focal plane. 



Acknowledgments 

We wish to thank Dennis Wright, Makoto Asai, and 
Daniel Brandt for their assistance and expertise in de- 
veloping Geant4 applications, and Gordon Bowden for 
useful discussions and critique. LSST project activities 
are supported in part by the National Science Founda- 
tion through Governing Cooperative Agreement 0809409 
managed by the Association of Universities for Research 
in Astronomy (AURA), and the Department of Energy 
under contract DE-AC02-76-SFO0515 with the SLAG 
National Accelerator Laboratory. Additional LSST fund- 
ing comes from private donations, grants to universities, 
and in-kind support from LSSTC Institutional Members. 



[1] Z. Ivezic et al., arXiv:0805.2366 (2008) 

[2] J. Singal, R. Schindler, C. Chang, C. Czodrowski, and C. 

Kim, Rev. Sci. lustrum., 81, 025101 (2009) 
[3] H. Howell, J. Wehrle, and J. Jostlein, Argonne Natl. Lab 

CP, 910505-191 (1991) 
[4] D.H. Davis, Journal of Applied Physics, 21, 1169 (1960) 



[5] S. Agostinelli et al.. Nuclear Instruments and Methods A, 
506, 250 (2003) 

[6] J. Allison et al., IEEE Transactions on Nuclear Science, 

53, 270 (2006) 
[7] P. Clausing, Ann Physik, 12, 961 (1932) 



